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ABSTRACT 

o 

' Higher derivative theory is one of the important models of quantum gravity, renor- 

(N ■ 

Qh| consider the 4 — e renormaUzation group for this theory, an approach which proved 

fruitful in 2 — e models. A consistent formulation in dimension n = 4 — e requires taking 
quantum effects of the topological term into account, hence we perform calculation 
which is more general than the ones done before. In the special n = 4 case we con- 
I firm a known result by Fradkin-Tseytlin and Avramidi-Barvinsky, while contributions 

0^ i from topological term do cancel. In the more general case of 4 — e renormalization 

^ , group equations there is an extensive ambiguity related to gauge-fixing dependence. 

As a result, physical interpretation of these equations is not universal unlike we treat 



, e as a small parameter. In the sector of essential couplings one can find a number of 

O ■ 

new fixed points, some of them have no analogs in the n = 4 case. 
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1 Introduction 

Quantizing gravitational field is an important aspect of quantum description of nature. Besides 
main efforts in the last decades were concentrated on the study of (super)string theories, recently 
there was an extensive development of a proper quantum gravity: both perturbative [1, 2, 3, 
4] and non-perturbative [5, 6, 7, 8]. It is well known that the perturbative approach meets 
serious difficulty which can be summarized as a conflict between renormalizability and unitarity 
[5]. Quantum gravity based on General Relativity is non-renormalizable, in particular higher 
derivative divergences emerge starting already from the one- loop level [1, 2]. In the presence of 
matter fields or at second loop [9] the divergences persist on-shell, indicating that no "magic" 
cancellations occur. The situation in supergravity extensions of General Relativity is somehow 
better, because higher derivative divergences emerge only at higher loops [10]. However, from the 
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general perspective there is no principal difference, for the corresponding version of supergravity 
extensions of General relativity are also non-rcnormalizable by power counting. 

An alternative way of quantizing gravity is to introduce some higher derivative terms into a 
classical action, treating them at the same footing as the lower-derivative (Einstein-Hilbert and 
cosmological) terms. For example, adding generic fomth order derivative terms, one modifies 
propagator and vertices in such a way that the new quantum theory is renornializable [3] (see 
also [11] for a more general formal proof of renormalizability) . This nice property enables one to 
establish the asymptotic freedom in the UV limit [12, 13, 4, 14] and explore the possible role of 
quantum gravity in the asymptotic behavior for GUT- like models [15]. Introducing new terms 
with derivatives higher than four leads to super-renormalizable theories of quantum gravity [16]. 

Unfortunately, the price for renormalizability is very high. The propagator of renormalizable 
quantum gravity contains unphysical poles corresponding to the states with negative kinetic 
energy or with a negative norm in the state space. These unphysical states are conventionally 
called massive ghosts [3] Despite the ghost masses have magnitude of the Planck order, they 
spoil unitarity of the S-matrix exactly in the range of energies where quantum gravity is supposed 
to apply^. Indeed the presence of massive unphysical ghosts in the tree-level propagator does not 
necessary imply that S-matrix is not unitary if one takes quantum corrections into account. It 
has been readily noticed [13, 23] that the quantum corrections to graviton propagator may, in 
principle, transform the real unphysical pole into a pair of complex poles. In this case the massive 
states should disappear in the asymptotic states and hence unitarity of the S-matrix gets restored. 
Unfortunately, until now there is no safe method to check whether this is really the effect of a 
resummation in the perturbative expansion. For example, despite the leading order of 1/N 
expansion provides corrections in appropriate form [23, 24], this approximation is not consistent 
in the case of quantum gravity [25] . Qualitatively the same is also true for other attempts in this 
direction and the final answer concerning the presence of ghosts in a full quantum theory remains 
unclear. In generic higher derivative gravitational theory (polynomial in derivatives) we meet 
complicated structure of the propagator, but with inevitable appearance of unphysical massive 
ghosts [16]. 

In order to clarify the situation with massive ghosts one needs to perform a non-perturbative 
treatment of the theory. In quantum theory, one of the known non-perturbative methods is the 
consideration of theory in a dimension 4— e , where e is a small parameter. This method has been 
developed by Wilson and Fisher for investigating critical phenomena [26] and later on applied 
to quantum field theory problems [27]. In quantum gravity the same idea has been consistently 
applied only in two-dimensional models [28, 29, 5, 30, 31], where it proved very fruitful. The 
main lesson of the two-dimensional case is the correspondence between 2 + e approach and 
numerical non-perturbative methods, such as dynamical triangulations [6]. One can accept this 

^General introduction to the fourth order derivative quantum gravity can be found in [17]. A detailed inves- 
tigation of propagator in higher derivative gravity models has been recently performed in [18] (see also further 
references therein). 

''There are interesting examples of higher derivative gravity theory with torsion which are free of ghosts [19, 20], 
but renormalizability can not be achieved in these models [21, 22]. 
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is a hint indicating the potential importance of the 4 — e consideration. In the present paper we 
shall consider the 4 — e approach for higher derivative qTiantum gravity. In particular, it looks 
interesting to explore the possibility of new fixed points in 4 — e theory. The important aspect is 
that the Gauss-Bonnet term, which is a part of higher derivative gravity action, is quite relevant 
for this purpose. The reason is that, for the theory in n = 4 — e dimensions, the Gauss-Bonnet 
term is not a topological term, in particular it contributes to many vertices and may be, in 
principle, affecting the /3-functions in the 4 — e case. 

The next question is whether the effect of the Gauss-Bonnet term has to vanish in n ^ 4 
limit. In fact, the opposite situation can not be ruled out beforehand [32]. The point is that 
the topological nature of this term is closely related to its diffeomorphism invariance. However, 
when theory is quantized (e.g., through the Faddeev- Popov procedure), the covariance is broken 
such that the vector space extends beyond physical degrees of freedom. After quantization, not 
only spin-2, but also spin-1 and spin-0 components of quantum metric become relevant, and 
the topological term creates new vertices of interaction between these components [32]. The 
quantum-gravitational loops may be, in principle, affected by the presence of topological term. 
Of course, this output does not look probable because, after including the topological term, the 
gauge-fixing condition should modify and eventually compensate the new vertices. But this is 
a sort of believes which are always good to verify and we will perform such verification in what 
follows. In a recent paper [33] we performed similar calculation for the particular case of Weyl 
quantum gravity and found that, in this case, the effect of Gauss-Bonnet term is relevant for 
4 — e renormalization group but vanish when n ^ 4. An extra benefit of the calculational 
scheme taking the Gauss-Bonnet term into account is the highest degree of safety for correctness 
of calculations. In case of Weyl theory we have confirmed the previous results [4, 34] for the n = 4 
renormalization group and also established conformal invariance of one-loop counterterms. Here 
we shall generalize this calculation for generic non-conformal theory and present a more complete 
investigation of 4 — e renormalization group, taking into account a gauge- fixing dependence and 
the corresponding ambiguity. 

In order to complete the picture and emphasize even more the relevance of the problem, let us 
remember that the Gauss-Bonnet-like term is an important ingredient of an effective low-energy 
action of the (super)string theory. The string effective action is composed by the Einstein-Hilbert 
term and an infinite set of higher derivative terms which emerge in the next order in a' (see, e.g., 
[35]). Therefore, models like effective low-energy quantum gravity [36, 37] which are supposed to 
be a kind of some universal low energy limit of string theories, may be sensitive to the presence 
of Gauss-Bonnet term. The uniqueness of the output in effective approaches depends on whether 
the n ^ 4 limit is universal or not. Hence in case we find vanishing contribution from the 
Gauss-Bonnet term in n — > 4 limit, this will be a certain confirmation of universality of quantum 
corrections. 

Thus, the main purpose of our work is to derive the one-loop corrections for higher-derivative 
quantum gravity with the Gauss-Bonnet term in n 7^ 4 and analyze the consequent renormal- 
ization group. The paper is organized as follows. In the next section we introduce notations 
and briefly discuss the theory. Section 3 contains brief summary of a Lagrangian quantization 
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of the theory and in particular the discussion of gauge-fixing conditions and related arbitrariness 
in counterterms. In section 4 we derive the trace of the coincidence limit of second Schwinger- 
DeWitt coefficient for an arbitrary dimension n. Furthermore, we consider the n ^ 4 limit and 
obtain the expression for the divergences. The cancellation of contributions from topological term 
and perfect fitting with the previous calculation of divergences in the higher derivative quantum 
gravity [14] confirm the correctness of our calculation. In section 5 we perform analytical and nu- 
merical investigation of 4 — e renormalization group equations for the higher derivative couplings 
and in section 6 consider the renormalization group equations for the Newton and cosmological 
constants. Finally, in section 7 we draw our conclusions. Many technical details and bulky for- 
mulas are collected in Appendices. Throughout the paper we use Euclidean signature in order to 
be consistent with the previous publications [4, 14]. 

2 General framework and notations 

The classical action of the theory has the form 

S = -f^""-^ j d"x^[ - + ^R' +TaR-^{R-2A)], (1) 

where // is some dimensional parameter, 

4 2 

= R^uapR'^ 7,R'OipR°''^ + 7 7T7 7^R^ 

n — 2 (n — l)(ra — 2) 

is the square of Weyl tensor, 

E = R^,,^pR>'-'^l^ - AR^pR^P + R^ 

is the integrand of the Gauss-Bonnet term, which is topological in ri = 4 and A, p, ^, r are 
independent parameters in the higher derivative sector of the action. Furthermore, □ = g^^V ^Vu 
and = I&ttG, where G is the Newton constant and A is the cosmological constant. 
An alternative form for the action (1) is 

S = -^(n-^) J (Tx^ Rl,^f^ + y Rl, +zR' + TnR - 1 (i? - 2A) | , (2) 

where wc denote Rf^^^^p = R^iyapR^'^"^ and i?^^ = RajjR"'^ ■ The parameters x, y, z are 
related to p, A, ^ as follows: 

_11 __2 4 __11 1 

''~2A"p' ^"~(n-2)A^p' ^ ~ "p^i^ A(n-l)(n-2) ' ^' 

The two forms of the action (1) and (2) are completely equivalent and we shall use one or another 
depending on convenience. 

At the classical level in the four dimensional space the role of the Gauss-Bonnet J £J-term 
is negligible. However, situation changes when theory is considered in ra / 4, where the j E- 
term becomes dynamical and hence must be taken into account. Using the language of Feynman 
diagrams, one can say that despite this term does not contribute to graviton propagator on flat 
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background for any n, it does contribute to vertices of interaction between different modes of 
metric (spin-2, spin-1 and spin-0). Hence the trivial nature of this term must be taken with 
proper caution, especially in n 7^ 4. Let us mention that, for n = 4 — e, there is a natural limit 
e = 1 imposed by the following circumstance. For n = 3 the expressions for E and coincide, 
and since the Wcyl tensor vanish at n = 3, the Gauss-Bonnet term also vanish. 

At quantum level, all terms in the action (1) are necessary for rcnormalizability. At the same 
time, the term has a key role. If we do not include this term, the quantum theory will 
have, in the graviton (traceless and transverse) sector of quantum metric the situation when the 
propagator has only second derivatives while there are four derivative vertices. Then the standard 
evaluation of the superficial degree of divergences for the diagrams of the theory (see, e.g. [17]) 
shows very bad renormalizability properties which are even worse than the ones for the quantum 
gravity based on General Relativity. On the other hand, since the term is included, the 

propagator has, in the momentum representation, the form 



where the negative sign of the second term indicates the ghost nature of the propagating spin-2 
state with mass m2- Other unphysical mode is possible in the spin-0 sector of metric, depending 
on the coefficient of the / i?^-term [3, 17, 18]. 

Without /2^-tcrm, the higher derivative sector of theory possesses (in n = 4 case) local 
conformal invariance and its rcnormalization has some special features. The main point is that 
conformal symmetry is violated by anomaly, therefore this symmetry is not supposed to hold in 
rcnormalization beyond one-loop level. We have recently discussed rcnormalization in conformal 
quantum gravity in [33] and hence will not consider it here. 

The relevance of the Gauss-Bonnet term at quantum level is also significant. Long ago has 
been noticed that the topological nature of this term is related to Bianchi identity [32]. As 
we already mentioned in the Introduction, it is in principle possible that the loop contributions 
depend on the presence of topological term. In particular, there is no reason to expect that the 
effect of Gauss-Bonnet term is negligible for n = 4 — e-dimension rcnormalization group, for this 
term is not topological at n 7^ 4 even at classical level. In previous paper [33] we have confirmed 
these considerations for the special case of Weyl quantum gravity and will now generalize these 
results for general higher derivative quantum gravity (1). 

3 Quantization and gauge fixing 

The general scheme of Lagrangian quantization for theory (1) has been formulated in [3] (see 
also [17] for more detailed exposition). The most useful method of practical calculations is 
through the background fields method and Schwingcr-DcWitt technique (see [38] for an extensive 
introduction including the higher derivative case). The application of these methods to the case 
of higher derivative gravity has some peculiarities which have been revealed in [4] (see also [39, 17] 
for more detailed pedagogical consideration). 




(4) 
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The background field method imphes special parametrization of the metric 

diJ-v 9f_iu ~ dfiv + hfj,i, . (5) 

In the r.h.s. of the last formula, g^i^ is a background metric and hfj_y is a quantum field (inte- 
gration variable in the path integral). The 1-loop effective action F^^^ can be written as [4] 

f(i)[5^^] = ^ In DetTi: - | In Det F"^ - i In DetHg^, (6) 

where 7i is bilinear (in the quantum fields) form of the action (1), taken together with the 
gauge-fixing term 

Sgf = /x"-^ j Xa Y^f" Xf3 ■ (7) 

Ti-gh is bilinear form of the action of the Faddeev- Popov ghosts, n is the renormalization parame- 
ter in dimensional regularization. The expression (6) includes also the contribution of the weight 
function y"^. In the case of higher derivative quantum gravity theory this term gives relevant 
contribution to the effective action, because V^/^ is a second order differential operator [4] . 

The general form of the gauge fixing conditions (here we restrict our attention to linear 
background gauges) has the form 

Yf,u = ^ (^g^^uO + tV^V^ - <5V^V^ + piR,,^ + P2Rgf^u) , (8) 

where = {a, (3, 7, 5, pi, P2) are arbitrary gauge-fixing parameters. The action of the Faddeev- 
Popov ghosts can be written as 

Sgh = /x"-' J (Px^ C-^ (Hghr^ a , (9) 

where 

Hgh = {-Hgh); = -S;a- V-V^ - 2/3V^V'^ . (10) 

Generally speaking, the counterterms may depend on the choice of gauge- fixing parameters, 
but this dependence is constrained by the on-shell gauge-fixing independence. Let us consider 
these restrictions, generalizing similar discussion of [40] for the general n case. 

We denote F(ai) the effective action corresponding to arbitrary values of gauge parameters 
and Trn = T{af^) calculated for special values of these parameters, a^^-* . Our purpose is to 
evaluate the gauge fixing dependence F(aj) or, cquivalcntly, the expression for the difference be- 
tween the two effective actions T{ai) — V„i- In general, this expression may be rather complicated 
(sec the details in [40]), but the part of r( o^'i) — F^ which is relevant for the renormalization 
group may be easily evaluated without special calculations. For this end we remember that the 
gauge dependence of counterterms has to disappear on the classical mass-shell (see, e.g., [41] for 
a general proof for gauge theories). Hence we can write 



T{ai) = Tm+ f (Px^ f^,{ai) , (11) 
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where fnv(pii) is some unknown function. The integration is taken over n-dimensional space, 

because our target is the rcnormahzation group in n = 4 — e dimensions. 

Since the object of interested is the divergent part of the effective action T{a,i), one can 
assume it is a local expression. Furthermore, both T„i and T{ai) have the same dimension as the 
classical equations of motion SS/Sg^i,. Therefore f^uioti) is a symmetric dimensionless tensor 
and the unique choice for it is 

fnA(^i) = 9ixv X f{ai) , (12) 
where f(ai) is a numerical quantity. Thus we arrive at the relation 

r(a,) = r„ + /(aOx fd^x^g^,^. (13) 

According to the last relation, the gauge-fixing dependence of divergent part (or, more general, 
the local part) of one-loop effective action is proportional to the trace of classical equations of 
motion. The corresponding traces for the relevant terms have the form 

flp^^ / "^'^ ^ = V5 { ^ R^axe - 2 ni?! , 

g,.^Jd^x^gR% = ^g{^Rl,-^nR^ , 

g,.-^ J d^x^{R- 2A) = ^|!i-^i?-nA|. (14) 

At this point it is better to separate the n = 4 and ra = 4 — e cases. For n = 4 we obtain the 
following formula for the gauge-fixing dependence [40]: 

r{ai) = Tm + f{ai) X J d^x^ { - ^ (i? - 4A) + | Di? } . (15) 

According to (15), the coefficients of and B? poles are gauge- fixing invariant. One of the 

consequences is that the corresponding renormalization group equations are universal, providing 
information about the UV limit of the theory. At the same time, (□i?)-type pole depends on 
the choice of the gauge-fixing condition. Taking different values of gauge-fixing parameters ai 
one can provide any desirable value of the function /(a^) and hence change the (□i?)-type 
counterterm. Thus, the corresponding parameter r in (1) is not essential. The immediate 
conclusion is that there is no much interest to calculate the renormalization of r, especially if the 
calculation is done for a particular gauge-fixing. Another situation takes place for the Einstein- 
Hilbert and cosmological terms. Renormalization of each of them is gauge-fixing dependent, but 
one can easily check that the dimensionless combination k^A is an essential coupling constant 
with the invariant renormalization relation. 

In the n = 4 — e case the situation is more complicated. Simple calculation using (14) yields 
the following generalization of (15): 

V{ai) = Tm + f{ai) X fi^-^ J d^x^ | [2x + ^ + 2(n - l)z\ (Di?) 
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+ + yi?^, - (16) 

An important consequence of the cq. (16) is that neither one of the parameters x,y, z,t, n, A 
is essential in the n = 4 — e case. However, gauge-fixing dependence is concentrated in a single 
numerical function /(a^) and therefore we can easily extract combinations of the couplings which 
are essential parameters. This result will be extensively used in sections 5 and 6, where we 
consider renormalization group equations for essential couplings. 

4 Bilinear expansion and derivation of divergences 

In order to apply the background field method, we need the bilinear expansion in /i^j, for the 
classical action (1). It proves useful to work with the equivalent form of the action (2), and 
therefore we expand this action as follows: 

- \ [v^i?-2VffA](2)|. (17) 

The expressions [^Rl,^^] , [y^i?^ J , [^R^] , [y/gR] [ A ] were derived in the 
previous work [33] and we will not reproduce them here^. Using these formulas and the expression 
for the gauge- fixing term (7), after performing some cumbersome commutations [33] one can find 
the bilinear form of the action 

[S + SGpf^ = h^'''n^.,aph''^. (18) 

The operator IrL = TL^y^ai} depends on gauge parameters, a, (3, 7, 5, pi and p2- For practical 
calculations these parameters have to be chosen in such a way that H assumes the most simple 
(minimal) form 

n = ka^ + o{v^), (i9) 

where is a non-degenerate c-number operator. The expressions for the gauge-fixing param- 
eters providing cancellation of the non-minimal four derivative terms g'^i/VaDV^ , ^a/jV^nVj^ , 
Vyi^Vj/VaV/j and ^i^/jV^nVa are the following: 

y + Az 2 2x-2z 

With this choice of the gauge parameters, the operator H takes the form (19) 

n = ka'^ + DP^VpVx + N^Vx - (VxZ^) + W , (21) 

where K, DP^, N^, and W arc local matrix expressions in the h^'^-space. The identity 
matrix in this space is defined as a symmetric tensor 

i = Sf^u,al3 = 7; {9lJ-a9u/3 + 9^13 gi^a) ■ (22) 



^AU but the first expansions can be also found in [12, 29, 14, 17]. 
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Before writing down the bulky formulas for the elements of (21), let us notice that the expressions 
for and are in fact irrelevant. The reason is that both of them arc proportional to the 
covariant derivatives of the curvature tensor. Therefore, due to the locality of divergences, these 
terms may contribute only to irrelevant gauge-fixing dependent J counterterm. Since we 

are not calculating this term here, for the sake of simplicity, in what follows we shall simply set 
both and to zero. 

The matrix coefficient K of the fourth derivative term has the following form: 



y + 4:X 



y + 4:z 



iK)^u,al3 = ^ <^|Uz/,a/3 + JJ^ — 9ixv Qaf) ■ (23) 



4(x - z) 



The expressions for Df^ and W have the form 



{D''^)i,i^,a/3 = -2x g^./3Ra''^^^ + 4x SPR^at^fi + (3a; + y) g^^R^ai^f} 

- {4x + 2y) 5P,g^p R\ - 2x g^pR^^gP^ + y g^,u5P,R^p 

- 2z 5apP^Rf,^ + 2z g^r^gP^Re^p + 2x g^^R^P^p + 2x5^pP^R^ 

+ (| ^ - ^) {^^t^y.c.P - g^lu gap g^^ - 2g^(3Sf^af^ + 2g^y5al3, 



+ ^^S^,,^pRP>^ - ^g^,gapRP\ (24) 
2>x + y y — bx y + 2x 

+ TT^R IX uRaarp "I" 2 R/xtR av(3 "I" 2 RfjhaRvP 

+ 2 ~ 4^2 ) (-^MO'^/S ~ 2 '^dvpRap) + zRfj,yRap 

- \ xRlxcT + yRl\ + ^^^ - ^i^- 2A) ( 5^1. , a/3 - ^ g^ugap ) 
3y 

+ -^di^f^^fJ-r^^a — Xg^uRaearRp " yg/iiyRarR^ (3 ■ (25) 

In the last two expressions we have used special condensed notations which help presenting them 
in a relatively compact way. All algebraic symmetries are implicit, including the symmetrization 
in the couples of indices (a/?) ^ (/UJ^), (o ^ /?) and (/i ^ i^) and in the couple (p ^ A) in 
the operator DP^. In order to obtain the complete formula explicitly, one has to restore all these 
symmetries. For example, 

RfipR^ aufS ^ 2 i.^fJ'P^^ ceufS "I" RapR^ p,Pv) 

restores the (a/?) {jiv) symmetry. Finally, let us notice that formulas (23), (24), (25) coincide 
with similar expressions of [14] (see the last reference therein) in the particular case l/p — 
after imposing the corresponding constraints on {x,y,z) in (3). 

The solution for the inverse matrix can be easily found in the form 
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where 

^ ^ y + 4z 

4x -Az + n{y + 4z) ' 

Let us remark that in dimensional regularization the divergences of the expressions In Det H 
and In Det {K~^ H) are the same, because an extra factor is a c-number operator. 

By straightforward algebra, we obtain the minimal operator in the standard form, useful for 
application of the Schwinger-DeWitt technique [4, 38] 

H = R-^ n = in^ + V^^VpVx + U . (26) 

It is important to notice that the new expressions 

VP>' = R-^DP^ and If = R-'^W (27) 

do not possess the symmetry {a(3) ^ (p-i^), while all other symmetries arc preserved. The 
expressions for the two matrices U = {U)^u,ai3 and = {V'^^J/iu ,ai3 are very bulky and we 
settle them in Appendix A. 

The algorithm for the one-loop divergences of a minimal fourth order operator (21) is the 
following [4] (see also [38] for an alternative, more systematic derivation): 

- In Det H = - / cTx^ tr lim a2{x' , x) = - ^ , (28) 

2 div (47r)^(n — 4) J x'-^x n — 4 



where 



x'-^x ' ' 90 ''^'^ 90 36 6 



- \Rp>y'^ + ^Ry''p+ ^y''py\+ yS''''^'^ ■ ^^^^ 

Here TZpx is the commutator of the covariant derivatives acting in the tensor h°^l^ space, 

7^pA = [Vp , Va] . (30) 

The particular traces of the expression (29) arc collected in the Appendix B. The result must be 
summed up with the contributions of ghosts In Dot iigh and the weight operator In Det Y"^ , 
according to the formula (6). The corresponding expressions are also given in Appendix B. Let us 
present just the overall result for the divergent part of the effective action, in terms of parameters 
A , p and ^ 

A\ = -/x"-^ j d^x^ I A(n) E + P2{n) + (iz{n) R" + ^ R 



The coefficients (/^-functions, as we shall see later on) are given by the expressions 



A(n) = 77^ 4°) + ^ + ^ , i = (l,2,3,4,5,6). (32) 



1 

(4<|=^ 
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For the sake of convenience let us now replace n = 4 — e, such that all coefficients , A) 

should be expressed in terms of e. The explicit form of the coefficients S^''^ corresponding to f3- 
functions for p, A and ^ are collected in Appendix C. Their general dependence on the couplings 
has the following structure. For /3i we have 

For P2 we have the expressions 



(0) _ .(0) C M) AO) A (1) _ (1) (1) 

2 — "2A °2B + ^2C ^ ' "2 — "2yl ? + "2^ ' 

5?'=5fle + 5?U^ + 5fix\ (34) 



The coefficients of /^s can be written as 



= 45 + 4°i 7 + 4S i , 4^^ = 4:i ^ + 4i^ a + 4g 



4^^ = 45e^+41?A+4gA^ (35) 



Furthermore, the coefficients 6^^ have the form 

4°^ = 4°ie+4°iA+42f > 4^^ = 4!ie^+42 A^ 4^) = o. (36) 

Finally, the coefficients of functions and /^g have the form 

4°^ = 45 e + 42 A , 4°^ = C e + 42 a^ , 4^-^^ = 4^'^^ = o . (37) 

The coefficients 4a BC depend only on e and are listed in Appendix C. 

After taking the limit n — 4 in coefficients /3j(n), we arrive at the counterterm, that is the 
negative of the n — ^ 4 coefficient for the pole term 

- (47r)2(n-4) 7 ^ \ ^ ^ I5" ^ 



^2 ^ ' 36/ " ' V12A 6 ^ 
3"~9A/"k2"^V72^~J^ 



The last expression does agree with the well known result of Avramidi and Barvinsky [14] ^. This 
coincidence is remarkable in several aspects. First of all, one can see that the effect of the Gauss- 
Bonnet term is not relevant for one-loop renormalization. This means that, in the framework 
of fourth derivative quantum gravity, the hypothesis of Capper and Kimber [32] concerning the 
relevance of the topological term on quantum level may be valid for finite corrections (e.g., this 
is the case for the 4 — e renormalization group) or for sub-leading divergences at higher loops 
but not for the leading logarithms, including the one-loop divergences. Since the theory under 
^Including the coefficient (5/36) , all coefficients agree with the previous calculation of [4]. 
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consideration involves all degrees of freedom of the quantum metric, it is likely that the same 
situation holds in any theory of quantum gravity. Second, the enormous cancellation of the 
p-dependcncc (sec Appendix C) in e ^ limit may be considered as a very strong test for the 
result (38) and [14]. Let us remark that Rcf. [14] corrected the result of previous pioneering 
calculations [12, 4], and therefore an extra verification docs not look unnecessary, especially in 
view of existing applications [42]. One has to notice that the difference between the results of [14] 
and [4] was only in i?^-term and that the work [14] included an additional test (derivation using 
C-regularization in a space of constant curvature) for some special combination of this and C^- 
term. That test could be inefficient only for the error proportional to a combination Ej^^ — ^E?. 
Our calculation represents a very strong test of this possibility and one may be certain that the 
result (38) is a correct one. Third, our calculations have been organized in such a way that we 
could separate the result for the conformal case | — > at every stage. Therefore, the coincidence 
of the final result with the one of [14] is an additional confirmation for our previous derivation of 
divergences in conformal quantum gravity [33], where wc met a perfect agreement with the result 
of [34] and partial agreement with the one in [4]. Finally, the cancellation of a p-dcpcndcnce 
in the e ^ limit provides a particular but rather strong test for the correctness of e 7^ 
coefficients presented in (33)-(37) and Appendix C. These coefficients will be applied in the next 
section for investigating the 4 — e renormalization group equations. 

5 Renormalization group equations for higher derivative terms 

Despite the p-dependcncc cancels out in the e ^ limit, the 4 — e renormalization group 
equations do not assume exactly the form of the known equations in four dimensions. The reason 
is that 4 — e renormalization group /3- functions are sensitive to 0(e)-corrections which depend 
on p. A standard derivation (see, e.g. [17]) yields the following renormalization group equations 
for the parameters of the higher derivative sector 

dp 1/9 2 

Tt = -'P + (^'P 

dt 

§ = -ei-M\ (39) 

where is the renormalization parameter of dimensional regularization, dt = dfi/fx and the 
/^-functions are given by formulas (32). In the limit e — ^ 0, we meet usual four dimensional 
renormalization group equations, which are exactly the ones obtained earlier in [14] 

/A \2dp 196 2 /A \2d'>^ ^ 2 

(47r)^§ = -lOX' + 5XC-^e. (40) 

It is easy to see that the n = 4 — e equations (39) are much more complicated than their 
4-dimensional cousins (40). In order to see this, let us notice that the first two equations (40) 
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do not depend on parameter ^. These two equations can be easily solved and their solutions 
replaced into the last equation, which can be also solved analytically [4, 14]. For convenience we 
present the known solution of all three equations and their analysis in Appendix D. However, 
the equations (39) do not admit any simple factorization and have to be solved simultaneously. 
Taking into account their complexity, there is no hope to achieve solution in the analytical form, 
and one has to apply numerical methods. 

In fact, the situation with 4 — e renormalization group equations is even more complex, 
because one has to account for the arbitrariness coming from the choice of a gauge-fixing con- 
dition. We have already considered this aspect of theory in section 3. Taking the gauge-fixing 
arbitrariness (16) into account, we arrive at the complete form of 4 — e renormalization group 
equations for the three parameters 

^ = -ep + epf{ai)+p^Pi, (41) 
^ = _eA + €A/(«,)-2A2/32, (42) 

f = -eC + eU{ai)-eP3, (43) 

where the /3-functions (32) correspond to minimal gauge-fixing condition. The remarkable feature 
of the gauge-fixing dependent terms is that they are proportional to the same quantity /(aj). 
Furthermore, the /3-functions (32) are homogeneous functions on the couplings. Taking all that 
into account, our main strategy in the investigation of the system (41) - (43) will be the following. 
We construct two combinations 9 and to of the effective charges (A, p, ^) such that the renormal- 
ization group cqTiations for these parameters are free from gauge-fixing ambiguity. The equations 
for the charges 9 and u; will be explored with the main purpose of establishing the UV stable 
fixed points. After that we shall consider the equation for the remaining effective charge, with 
the invariant combinations and lo at the fixed point. In this way we shall learn the asymptotic 
UV behaviour corresponding to a given fixed point. 

The most natural choice for independent eflFective charge is of course A, for it defines the 
interaction between gravitons. Therefore we define the invariant charges as ratios between other 
parameters and A 

= -, ^ = -^- (44) 

The coefficient in the second expression provides correspondence with the notations of [4, 14]. It 
is straightforward to see that the renormalization group equations for 9 and co are independent 
on function /(ctj) and have the following universal form: 

^ = -2ep2-pi, ^ = -2a;^2-3/33. (45) 



Here T{t) is a new parameter defined by 

m 

(4vr) 



dr = TTzk , (46) 
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where A(t) is a solution of (42). 

Equation (46) can be viewed as a most important relation of the whole 4 — e approach. We 
know that t ^ oo corresponds to the high energy limit ^ oo in the MS renormalization 
scheme. However, due to the gauge-fixing dependence of the renormalization group equation 
for A(t), it is not clear whether this high energy limit corresponds to a certain limit in the new 
variable r. Perhaps even more important aspect of the same problem is the possible arbitrariness 
in the asymptotic behaviour of the coupling A. This arbitrariness may put under question the 
asymptotic freedom of the theory, spoiling the physical sense of the renormalization group. We 
shall consider the asymptotic behaviour of A later on and for a while, when looking for the new 
fixed points for uj and 0, will identify the limit r — oo as the UV one. 



Fixed Point 


6 




Stability 


1 


0.33516 


-5.38892 


Saddle 


2 


4.61183 


-1.60198 


UV-Unstable 


3 


-4.31710 


-1.47066 


UV-Unstable 


4 


-4.44192 


-0.15162 


Saddle 


5 


4.80565 


-0.06229 


Saddle 


6 


0.33782 


-0.00283 


UV-Stable 


7 


-3.911(32 


0.03123 


ITV-Uustal)lc 


8 


-4.11072 


0.07230 


Saddle 



Table 1. The list of the fixed points and their stability with respect to small 
perturbations for the case e = 0.1. 

According to the outline described above, consider equations (45). For any given value of 
e, this is a system of two nonlinear and rather complicated equations. There are no chances to 
obtain the general solution analytically, so the unique visible possibility is to fix the value of e, 
according to the standard practice [27]. We shall try different values of e, starting from e = 0.1. 
In this case the equations for the fixed points 

2^/52 + /3i = 0, 2u;/32 + 3/33 = 0, (47) 

have eight distinct solutions collected in the Table 1. We investigated their stability under small 
perturbations of both charges Ok Ok + SiOk, tOk ^ tOk + ^^k^ where the values Ok, iOk, 
k = 1,2, ..,8 correspond to different fixed points. It turns out that for e = 0.1 there are all 
kinds of the fixed points: UV stable, UV-unstable in all directions and saddle points. 

In order to compare with the standard e = case (see Appendix D), we present the 
corresponding fixed points in the Table 2. In this case there are only two fixed points, what 
means that six new fixed points have emerged due to the procedure e = 0.1. The renormalization 
group trajectories which result from numerical investigation are shown at Figure 1 for the simplest 
e = case. Each arrow indicates the slope of integral curves, dO/du at a given point of the phase 
diagram. 
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Figure 1: Diagram for e = (four dimensions). Here y stands for 9 and x for 



Fixed Point 


e 


u 


Stability 


1 


0.32719 


-0.0228(3 


UV-Slal)le 


2 


0.32749 


-5.4671 


Saddle 



Table 2. Fixed points for 9 and uj and their stability for the simple case e = 
The equations admit also complete analytical investigation presented in Appendix D. 

For e = 0.1, the renormalization group flow is illustrated at Figure 2. Notice that the points 
1 and 6 are shown in detail in the right diagram, which is very similar to the diagram at Figure 
1. Indeed the trajectories in a limited region containing points 1 and 6 are practically the same 
trajectories carried out for e = 0. Thus, the procedure e = 0.1 keeps the dynamics unchanged 
in some limited region of the plane^, producing just a small displacement of the fixed points (in 
this sense, points 1 and 6 are direct descendants from the standard 4D-fixed points). However, 
in large scale, dynamics is totally modified by the appearance of six extra fixed points. 

The situation for other small positive values of e is qualitatively similar to the one in the 
e = 0.1 case. For example, for e = 0.01 we meet 10 fixed points, but only one of them is 
UV-stable. We will not present the details of this case here. 

The results for the negative e are quite different. The result of numerical integration for the 
cases e = —0.1 and e = —0.01 are shown at Figure 3. The general structure is similar to the 
solution in four dimensions, and no extra fixed point emerges by consequence of e 7^ 0, only small 
displacements of these points take place. 

Now we are in a position to explore the issue of asymptotic freedom in 4 — e framework. For 

this end we need to investigate the behavior of the effective charge A in the UV limit t — > 00. 

Also, this defines the physical sense of parameter r(t) and its possible relation to the change of 

^The meaning of this statement is somewhat restrictive: of course one has to consider only trajectories confined 
in that region. 
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Figure 2: Numerical integration for e = 0.1 in the whole plane (a;, ^) in the left diagram. The 
right diagram shows a more detailed description of the trajectories between point 1 (saddle) and 
point 6, which is UV-stable. 

a physical energy scale. Hence, this consideration should clarify the physical sense of equations 
(45) and their main characteristics, that are the corresponding fixed points. 

Starting from the equation for A(f) in (39), one can find the analytical form of A(t) in the 
vicinity of a fixed point (wQ) ^o)- For this end, the expression (i^ must be rewritten in terms of 
and which should be further replaced by cuq a-i^d Qq. After performing this, independent on 
the values of and we arrive at the equation 

^ = aA-62A2, (48) 

where 

a = -e + e /(a^) and }? = 2 , ^o) • (49) 

Let us remark that the parameter h in the last equation depends only on the values (wqj^o)) 
while the parameter a depends also on the choice of a gauge-fixing condition and therefore can 
be made arbitrary. 

The solution of (48) is straightforward 

Starting from this relation, we integrate (46) and arrive at the explicit form of r 

r(i) = ^ In [a + 6A0 (e«* - 1)] + C , (51) 

where C is irrelevant integration constant. 

It is easy to see that, formally, the asymptotic freedom in the theory depends on the sign of 
the quantity a in equation (49). In case a > we have A(t) — in the UV t — 00, and also 
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Figure 3: Numerical integration for e = —0.1 (left) and e = —0.01 (right), where the two fixed 
points in both cases are easily recognizable. 

r — oo, so everything is consistent with the e = case. On the opposite, for a < we meet the 
asymptotic behaviour A(i) — Aq 7^ in the UV limit i — ^ 00, also r — ^ tq 7^ 00. 
It is important that the explicit form of the function /(a^) derived in [40] 



3 A 



2 2 



97A5/2 



f{ai = a,7) = - <^ - In TT^^ ^ + a - \ ) (52) 



2a3/2 (a-6A7) 



does not produce any restrictions for the value of f{ai), which can be modified arbitrarily by 
choosing one or another gauge- fixing condition. According to equation (49), this means that 
one can also change the sign of the quantity a, and thus the asymptotic freedom in the 4 — e 
theory is not a physical phenomenon, but an artificial occurrence depending on the choice of the 
gauge-fixing condition. 

Is it true that 4 — e renormalization group does not have physical sense for the case of higher 
derivative quantum gravity? In fact, this is a problem of interpretation. Let us remind that the 
original 4 — e approach [27] treats e as a small parameter of the non-perturbative expansion. If 
we take this position and consider all terms proportional to e as small by deGnition, then the 
asymptotic behavior for A is close to n = 4 renormalization group equation (79). In this case 
the status of the new fixed points becomes clear and we conclude that the stable fixed points 
in all cases are similar to the one in n = 4 renormalization group. The same concerns also 
the renormalization group flows in the vicinity of the stable fixed points. At the same time, the 
existence of the numerous new unstable and saddle fixed points for e > indicates the possibility 
of a rich non-perturbative structure of the theory. 
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6 Renormalization Group in the Einstein-Hilbert sector 



Consider the renormalization group equations in the Einstein-Hilbert sector of the theory, includ- 
ing the cosmological constant. The renormalization group equations for the parameters in the 
low energy sector of theory are defined by the expressions for divergences (31), but must also 
account for the gauge-fixing dependence (16). The form of these equations is as follows: 

^ = („ _ 4) «2 + «2^4 + ^kV(«^) (53) 
for the Newtonian constant and 

^ = A/34 + ^ A/35 + ^ - A f{ai) (54) 

for the cosmological constant. 

The presence of the function f{ai) in both equations indicates a strong gauge-fixing depen- 
dence. Indeed, there is nothing new in this dependence, as already known from [4] (see also [40]). 
The standard solution of this problem is to look for the gauge-fixing invariant combination of the 
two /3-functions. It is assumed that the corresponding effective charge is an essential coupling 
constant [5]. Indeed, in the n = 4 case such combination is the dimcnsionlcss product k^A of 
the two effective charges k and A. The qTicstion is whether the essential effective charge can 
be found for the case of 4 — e renormalization group equation. 

By direct computation one can easily obtain the unique combination 

' 2 

of K and A which is explicitly independent on gauge-fixing parameters. As one could expect, 
this is exactly the combination of k and A which is dimensionless for any value of n. Of 
course, for n = 4 we come back to the standard combination of parameters 7(4) = k^A. 
Using (53) and (54) we arrive at the renormalization group equation for 7 



It ^ ^ 



If N 
n - 4 + - [2{N + l)/34 + m + ^/36 



(55) 



The last equation looks like the one for the essential coupling, but in fact there is still a difficulty 
in its interpretation. The problem is that the expressions /34,5,6 can be written in terms of 
essential couplings, u and 9, but there will be some extra factors of \{t) in (55) such that it can 
not be absorbed in a redefinition of dt. Besides, the gauge dependent quantity k^K (which is 
dimensionless only in four dimensions) also appears in the coefficient of (3q. The consequence is 
that, without imposing the smallness of e, the (4 — e)-rcnormalization group for k^A^ is gauge 
dependent, exactly as in the case of the higher derivative couplings. On the other hand, if we 
remember that all our consideration is called to model the non-perturbativc e-expansion and 
define e as a small parameter, the renormalization group flow near the UV-stable fixed point 
u;o, ^0 will not have qualitative difference with the e = case investigated in [4, 14]. The main 
feature of the UV asymptotic behaviour of 7 is strong dependence on the initial data for the 
couplings ^ and 7 . 
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The most interesting aspect of the renormahzation group equation for the cosmological con- 
stant would be to study the low-energy limit. This investigation could provide better under- 
standing of the possible role of quantum gravity for the cosmological constant problem [46, 37], 
including the remnant quantum effects of the particles with the mass of the Planck order of 
magnitude [47] and also for establishing the assumed universality property of quantum gravity 
[36]. However, the framework presented above is not appropriate for this purpose, because in 
the low-energy region we expect to meet a decoupling of the massive (ghost) mode (4) and the 
renormalization group equation should essentially modify. In order to observe decoupling one 
has to perform calculations and find /3-functions in the physical mass-dependent renormaliza- 
tion scheme. Similar program has been realized recently for the massive matter fields on curved 
background [48]. The practical calculations in higher derivative quantum gravity, despite they 
are looking rather difficult from the technical point of view, represent a serious challenge for the 
future work. 

7 Conclusions 

We performed a complicated calculation of the counterterms in the fourth derivative quantum 
gravity. It is the first time that the quantum effects of the Gauss-Bonnet term have been taken 
into account. In the n = 4 — e case the effects of this term are shown to be non-trivial, in 
accordance with the prediction by Capper and Kimber [32] concerning an important role of the 
topological term in quantum gravity. At the same time, for n = 4 the effects of the topological 
term cancel, in this limit we recover the known result [4, 14]. This coincidence represents very 
efficient verification of previous derivations of the n = 4 renormalization group equations. 

The renormalization group equations for the parameters of theory in the n = 4 — e are essen- 
tially more complicated than in the n = 4 case, in particular, the equations for different couplings 
do not separate. Moreover these equations manifest much stronger gauge-fixing dependence than 
in the standard n = 4 case. We separated the universal corner of theory, which is characterized 
by a number of new fixed points for e > 0. These fixed points may be related to the rich non- 
perturbative structure of the theory in the e-expansion. Out of the small-e approximation the 
equation for the most important coupling A (the parameter of the loop expansion of the theory) , 
manifests a strong gauge-fixing dependence which does not take place for the usual n = 4 renor- 
malization group. Prom a purely formal point of view, the physical interpretation of the n = 4 — e 
renormalization group looks unclear in this case. However the situation changes if we treat e as 
a small parameter. Then the theory remains asymptotically free for e ^ Q and moreover, as we 
already mentioned, has a number of new fixed points for e > 0. 
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Appendix A 

In this Appendix we collect bulky elements of the operator (26). The expression for U = {U)^v,al3 
has the following form: 

{U)^v,al3 = y _|_ 4^ I ^ gupRjj.pXcjRa'^^'^ H ^-^R^an'^Rxi^fip -\ ^^R'^a/j.^RuPpX 

+ (^RuaR'^au/S + RaaR" npi^j + ^~~2 ^t^'^^^l^ ~ -^daliRnecTRv^'^^ 

1 2 / r o \ ( ^^R ^ A D dX r> 

— ['>fJ.u,al3 - i'l gp-vQafi) + I — ^ I guf3^a^l + —gu^R p.RaX 

"I" ^ gixv ^ J^aajiT "I" I ~^ 4/^^ ) Rl^Oivji H ^ -'^ ;u vJ^aajS 

~ ^gapRti'^RuT — xi^i gp^uRaearRtf"^ — ^2 g^uRa^ R/3t + zR^i^Rafj 



+ ( ^ ~ ^ ) dapRiiv + [-^^Z- zR^i ) gf,uRal3 



where 



ki 


= gi^p g''^ Rpux , 


k2 = 5fj_u^ocl3g''^ , 


k3 = 


~ g^ RfxavP 1 


k4 


X p 

— "up, ^iia ) 


ks = ^vpf^ g^a , 


k6 = 


^IMv, R^ ) 


ky 


= \{5i'R^\.p, 


~\~ R ^ fii/a ) ) 


k8 = 


9upS\;r'^^), 


kg 


= gvP R{a''^ij.) , 


klO = 2 {^ajif^ Rfiv 




f^Ra/s), 


kii 


~ g^ivRcx^ p 1 


kl2 = gapRp.'^^u , 


ki3 ■■ 


= gp.v g^^Rafs 


ki4 


= gaf3 g^^R/j-u , 


kl5 = gp,u Rp , 


ki6 = 


= gaf3S^Ri, 


ki7 


= g^u , 


kl8 = ga/3 Sfj,t,P^ , 


ki9 = 


= gij.v gap g''^ , 


k20 


= g/xi^ gap R"^ 









(56) 



= xR%,, + yRl, + zR^-^{R- 2A) , (57) 



and the coefficients ^1,2,3 arc given by 

2x + 3y + l0z 7 xy + 9y'^ - 4x z + 28y z 

1 S ' ^ 45] 

n3 = ^ — , with E = Ax-Az + n{y + 4:z). (58) 

Zj 

For the matrix t^'''^ = (V^''^)//i/,a^ we have the expression 

. 20 

where the following condensed notations have been used: 



(60) 
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and 



zR 1 

bi = -2x , b2 = ^ , 63 = 3x + y , 64 = 2a:; 

1 y 

h = 2k?~ bQ = x + -, 67 = -4a; , bs = -4a; - 2y , 

69 = — 2x , 610 = —2z , 611 = 4x f^a , 612 = a; , 

^'13 = -y ^^3 , bi4 = z, 615 = 2y , = - , 

^ Oi-2J] zi? 1 

617 = 2zRn3 2^^' ^18 = ^-^, 619 = -617, 
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-y^3. (61) 



The last expressions have a non-symmetric form, and the rule of symmetrization is the one 
described in the section 4. 

Appendix B. Particular results for necessary traces 

In this Appendix we collect some intermediate formulas, which are necessary for the computation 
of the one-loop contributions. The calculations has been verified by using the package MathTensor 
[45] driven by the software Mathematica [44]. 
The trace of the operator U is given by 

ticU = 5^^'"^ U^^^ap = A, Rl^^p + A2 Rl, + A3 ii^ + ^4 4 + ^5 4 ' (62) 



where 



and 



. Pkl +Pk2n+ Pk?, +PkAn^ , ^ n q a k 

Ah = — 7z r-TZ : : — Ti fc=l,2, 3, 4, 5 

2 (4a; + y) (4x-4z + ny + 4nz) ' ' > ' > 



pii = S{8x'^ - 20xz - 3yz) , pu = 2 {l2x'^ + ISxy + 3y'^ + 4:0xz + 12yz) , 

Pi3 = -X {Ax -by -2Az) , = -x {y + Az) , 

P21 = SiUx"^ - Sy"^ - 8xz - 12yz - Az"^), P22 = 2{2Axy + 5y'^ + A8xz + 12yz + 16z'^), 

P23 = -y {Ax-5y-2Az) , P24 = -y {y + Az) , 

P3i = 8 (2x^ + xy- Axz - Ayz - lOz^) , ^32 = 2 {2xy + + 20a;z + 7yz) , 

P33 = -2; (4a;- 5y- 242;) , ^»34 = -^^ (y + 4^;) , 

P41 = 8{x + y + Sz), p42 = -Aiy + z + 3x), p4^ = Ax - 2y - I2z, pu = y + Az, 

P51 = , P52 = A{y + Az), p53 = -2 (Ax + y) , p54 = -2{y + Az) . 

The expression for tr (VyR) can be presented as 

tr {V^R) = BiR^ + ^R, (63) 

where 

^ _ ki + li2n + ks IT? + Iji + li5 -12 
~ 2{Ax + y){Ax-Az + ny + Anz)' ~ ' 
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and 



In = -16(4x-j/-2z)(2x + 7/ + 2z), 

/i2 = -4 (16x2 + 20xy + 5y^ - Axz + 6yz + 12z^) , 

/i3 = -16x^ -8xy-10y^ -72xz-56yz-8z'^, 

li4 = -4xy + 2y'^ - 8xz + 6yz - 16z'^, I15 = 2z (y + 4z) ; 

hi = , h2 = -4(2a; + 2z + y), hs = 4(3x + z + y) ^ 

hi = y - 4a; + 82; , h^ = -y - ; 

The computation of tr {Vf^Rpx) yields 

tr {VP^Rpx) = Ce R% + i?' + C3 4 



where 



^ _ Pii +Pi2n + pi3 + ■ - fi 7 

' ~ {Ax + y) {Ax-4.z + n{y + Az))' ' " ^' ^ ' 



C3 = 



(;?- - 3 + 2) ( i .c + ,y + 2 ,y + i + i n z) 
2 {Ax + y) {Ax-4z + n {y + Az)) 



with 



P6I = -8 {Ax-y-2z) {2x + y + 2z) , 

P62 = -2 (4a;2 + 14xy + 3y2 + 20xz + 8yz + 1622) , 

P63 = Sx^ + 2xy - 3j/2 - 16x2; - 16j/z, p64 = (2x + y) (y + 4^;) , 

P71 = —4 (6x + 3xy + y — 12 xz — yz — 2z), 

P72 = -2 {8x'^ + 3xy + y'^ + 10xz + 6yz + 2z'^) , 

P73 = —Axy — \2xz — yz — 8z^ p74 = z{y + 4z). 

The results for the traces ^ tr (F^^V/) + ^ tr (YxpV^^) can be written as 
1 tr {V^Vp + 1 tr = D, RI^^^ + Z?2 4. + D^R^ + D^^ + D, 

where 



96(4x + yf (4x - 4z + n (y + 4z)) 



^ _ gg + ft2 n + g.3 + + 9,5 + g.6 + m 
96(4x + y) {ix-Az + n (y + 42)) 

with ah non-zero given by 

qii = 1536 X (2^2 — 6x2; — y 2) , 

qi2 = 192 (26x^ + 18x^y + 3xy^ + 6x'^ z-2xyz-y^ z) , 

qi3 = A8 {I8x^ + 36x'^y + 12xy'^ + y^ + 78x'^ z + 28xyz + 2y'^ z) , 
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qu = 24(3x + y)2 (y + 4z) , 

^21 = -2bQ{%x^ + 2x'^y^-Axy^-y'^ + Ux^z + ll2x'^yz + 2{)xy^z 

- imx'^ -SQxyz^ -IQy'^ + ?>2xz^ -^z"^) , 

922 = 128 (I20x^ + 84x^y + 27x'^y'^ + 6xy^ + / - 224a;^z - USx'^yz - hAxy'^z 

- 6y^z - ISGx^z^ - USxyz'^ - SSy'^z'^ + 80xz^ - \2yz^ - 32z^) , 

923 = 16 (448x^ + IQAx^y + 360x^^2 + 64xy3 + 15y^ + 1024x3^ + 576x^^2 

+ Sxy^z + My^z - IMAx^z^ - 816xyz^ + llGy^z^ + 6Axz^ + lUyz^ + 128^) , 

924 = 8 (64x^ + 448x^y + SUx'^y'^ + 104xy^ + + 1664x^z + UOSx'^yz 
+ 560xy'^z + 28y^^ + 320x^z^ + bUxyz"^ + Sy^z^ - 512xz^ - 128yz^) , 
= 8 (y + 4z) (32x^ + 56x^y + 28xy^ + 3y^ + 192x^z + 112x^2 

+ 32x^2 + l2?/2^ + 8yz^) , 926 = 32x^(^ + 4^)2, 

931 = -128 {56x* + 32x''^y + ISx^y^ + 6x/ + / - 144x^z - 48x'^yz 

- 2y^z - 24xy^z + 240x^2^ - Gy^z^ - 80xz^ - IQyz^ - s/) , 

^32 = -64 (68x^ + 96x^y + 36x^^2 + 13xy^ + 3/ + 72x^2 - A8x^yz + lOxy'^z 

+ 13y^z - 216x^^2 - 16x^2^ + 26y'^z'^ + 168xz^ + 64^2^ + 36/) , 

^33 = -8 (32x^ + 304x^y + 204x^^2 + 44xy^ + y^ + 1088x^2 + 784x2yz + 280x^^2 

+ 4y^z - 1088x2 - SOxyz^ - 144^^^^ _^ 192^;/ - 336^2;^ - 160z^) , 

^34 = -4 (32x^y + lOOx^y^ + 32xy^ + 7y^ + 192x^z + 640x^^2; + 208xy2^ 

+ 1392x^^2 + QSy^z + A48xyz^ + 216^^2:2 _ 544x/ + 64yz^ - 80^) , 

qr35 = 4 (-4xV - 4x?/^ + - 64x2y2; - 24x2/^2 + ey^z - 144x^^2 

- ISy^z^ - 8xyz^ - 136yz^ - BOz'^) , 

936 = -4z (y + 4z) (4xy - 2y2 + 8x2; - 9yz + 42;2) , 937 = 4^2(^ + 4^)2, 

941 = -128 (8x^ + 4x2y + 6xy2 + y^ - 80x2^ - 8xyz + 40x^2 + ^yz^) , 

942 = -32 (I6x^ - 8xy2 - y^ + 120x2^ - 2y2z - 160x^2 - 32yz2 _ 8^3) ^ 

943 = 16 (iGx^ - 8x2y + 16xy2 + y^ - 176x22: + 24xy2 - 14y22 
+ 80x^2 - 48y22 - IQ2^>) , 

944 = 8 (lax-*^ - 12xy2 + y3 + 120x2^ - IGxyz + 16y22 - 160x22 - 402^) , 

945 = 4 (I6x2y - 4xy2 + y^ + 16x22 - 40xy2 + 18y22 + 88y22 + 802^) , 

946 = -4 (y + 42) (-2xy + y2 + 5y2 - 422) , 947 = -42(y + 42)2 , 
951 = , 952 = 16 (6x + y - 22) (2x + y + 22) , 

953 = -8 (24x2 + 12xy + y2 - 4y2 - 822) , 

^54 = -8 (6x2 _^ + y2 + 20x2 - 1022) , 

955 = 2 (24x2 + 12xy - y2 - 20y2 - 4022) , 

956 = (8x + y - 42) (y + 42) , 957 = (y + 42)2 . 
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Finally, let us present details of the contributions from Faddeev-Popov ghosts and weight 
operator (sometimes called third ghost). The action of the ghosts has the form 



Sgh = J <fx^ {ngh)l a , (67) 
where the operator Hgh is given by the expression 

Ugh = ingh)1 = ■^^^'^'P- 

Here V^^'p is generator of the diffeomorphism transformations for the metric 5g^^ = V^" 'p Sx^ 

T^liu ,/3 = -9fj,l3 - guf3 . 

By direct computation one derives, assuming pi = p2 = in (8) and keeping other gauge-fixing 
parameters arbitrary (10) 

Hgh = {ngh)l = -5;: □ - V^V^ - 2^V^V'^ . (68) 

In order to evaluate the expression —i In DetHgh, let us rewrite (68) in the form suitable for 
the generalized Schwinger-DeWitt method [38] 

T^gh = -[d]:a-a V^V + P;] , (69) 

where a = —(1 + 2(3) and P^i, = R^i,. The n dimensional analog of the known algorithm 
[4, 38] for the non-minimal Abelian vector operator (69) has the form 



2 



where P = and 

a _ 1 + 2/? 
^~ l-a~ 2 + 2(3 
Taking into account P^v = Rnu, we find 



(47r)2(n-4)7 "^^l 90 



-i In DetT^g^ldiv = jjz'^ ~ I ^^^Vd \ nn ^I'^c.p 



/I ,n 4, 90-n\ n /I ,n 1, n + 12\ ol 
+ (3"^ +3*'+-90-)^"''+(6"* +3^ + ^)^T 

The contribution of the weight operator (8) is calculated with the choice (20). The divergent 
part of Tr Iny can be computed directly from (70). As usual [4, 38], in the case P^i^ = —Rfxu, 
the gauge dependence is canceled (this cancellation is related to the gauge invariance of the 
electromagnetic field, whose contribution is given by exactly the same operator): 

-2 = (4.)^(n-4) J {^^-/^ + ^^^/^ + • (^2) 
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Appendix C. Coefficients of functions in 4 — e dimension 



This Appendix contains the coefficients 6 of the /3-functions (32), most of them were not included 
into the main text. The total expressions for the /^-functions is as follows (32) 

s P P P P ^ 

+ 45^+42^+4gf p3) 

The coefficients 5^*]^ ^ ^ have the form 

^(0) _ 37632 - 92704 e + 81860 - 31350 + 5533 - 446 + 15 
~ 2880 (e - 3) (e - 1)^ 



^(0) ^ e (e - 6) (6 - 2) 

48 (e - if 
^(1) _ e (e + 4) 



2 



12 (e- 3)2' 



e (e - 2) (92 - 176e + 57 - lOe^ + e^) 



48 (e - 3) (e - 1)^ 

(2) ^ (6-2) (6-1) 

96(6-3)^ ' 
^(2) ^ (6-5) (6-4) (6-2) 

48(6-3)3 



= ^ "'7or"2r ':r"^^^ (74) 



6^ (e - 5) (6 - 2f (6^ - 12 6 + 29) 
48(e-3)2(e- 1) 

For the expressions 62]^ b c ^ have 

6 (e2 + 4e-4) 



32(6-3)2 (e-1) ' 

(6 - 2) (-9576 + 21180 e - 15410 + 4410 6^ - 399 + 5 6^) 



^(0) 



^(0) = 

2^ 960 (e - 3) (6 - 1)2 

^(0) ^ 6(6-6) (6-2)2 

2^ 48 (6 - 1)2 ' 

^(1) ^ 6(6-2) (26 + 1) 

2^ 24(6-3)2(6-1)' 



6 (6-7) (6-2)2 ^^2_4g + iQ) 



2^ 48 (6 - 3) (6 - 1)2 



25 



(2) ^ (6-2) (e-1) 
96(6-3)4 ' 
^(2) ^ 6^6-2)^ (^2-86 + 19) 
48 (6 - 3)3 (6-1) 

^(2) ^ 6(6-2)=^ (_162 + 95 6 - 18 6^+63) 

48(6-3)2(6-1)2 ■ ^ ' 

The coefficients (Jg*]^ ^ q can be written as 



5. 



(0) _ - 1200 + 8808 e - 12480 + 5980 - 1290 e"^ + 127 - 5 

'3A ~ ;;:;r7 ^ 



960(6 - 3)2 
(6 - 2)2 (180 + 28 6 - 16 6^ + 6^) 
48 (6-3) (6-1) ' 



3 



^(0) ^ (6-6) (6-5) (6-4) (6-2) 

3^ 96(6-1)2 

^(1) ^ 6 (6 - 2) (6 - 1) (62 - 20) 

3^ 96(6-3)3 

^(1) _ 6 (6-2)2 (_io^246- 9e2 + e3) 

3^ ~ 24(6-3)2(6-1) 

(1) ^ 6 (( - 6) (f - 5) (f - 2)3 
48(e_3)(e_i) ' 

^(2) ^ 62 (6-5) (6-2)2(6-1)2 

3^ 192(6-3)^ 

_ 6^ (6-5) (6-2)2 

3^ " 48(6-3)4 



4c = ^ :^3^ ' (76) 



62 (6-5) (6-2)3 (6^-86 + 19) 

96(6-3)3 

Tfie coefficients ^4*]^ ^ q are the following: 

(6- 2) (72 - 406 -462 + 6^) 

192 (6 - 3) 



^(0) 



^(0) 



(6 - 2)2 (- 156 + 144 6 - 27 62 + 63) 
96 (6 - 3) (6 - 1) 



(5' 



;(0) ^ (6-6) (6-5) (6-4) (6-2) 

'^^ 96(6-1)2 



3 



6 (6 - 6) (6 - 2)2 (6 - 1) 



4A — „^3 



192 (6 - 3)' 

.(1) _ 6(6-6) (6-5) (6-2)3 
- 96(6-3) (6-1) ' 



The remaining coefficients corresponding to 5^^\ ^ q and 5q\ ^ q have the form 

e - 4) (6 - 2) .0) _ (6-4) (6-2) (6^-86 + 14) 



(0) _ V- -/ V- . _ 

4(6-3)2 ' 2(6 -3) (6-1) 
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^(0) ^ (e-6) (6-4) {e-2f ^(p) _ (e - 6) (e - 5) (6-4) (e - 2f 

768(6- 3)2 ' 6B 384(6- 1)2 ' ^ ' 

Appendix D. Renormalization group equations at n — A 

For convenience we present the solution [4, 14] of the n = 4 renormahzation group equations (40). 
The first two equations may be solved immediately. 



l + 62poi' 45 (47r)2' 

In the UV limit t — oo the solutions (79) manifest the asymptotic freedom for the two parameters. 
Indeed, the term "asymptotic freedom" is proper only for A, which is a parameter of the loop 
expansion of theory (1). The positivity in the theory is consistent with A > and therefore, 
according to (79), we need Aq > too. In contrast, is the coefficient of the topological 

Gauss-Bonnet term, which does not correspond to any gravitons interaction in the n = 4 case. 
The third equation (40) can be rewritten in terms of a new variable w{t) = ^{t)/\{t) 

dw C , , ,x C 



I* + ^ J = ^ - + ^) = ^ - - ^2) , (80) 



where 



^ = W?' 133' ^=133Tl8- -V2=o±\/--^- (81) 



Let us notice that w\ k l/k > 0. The solution of Eq. (80) can be easily found and 

written in terms of the original variables, notations (81) and wq = w{0), as 

w = ^^^^^ ^ ^ f \ m 

I- X \Wo-W2 J ^ ' 

and 

m = — ^ 5 0.517. 

The UV behavior of w{t) depends on wq. 

i) For Wq < W2 the UV limit is w ^ -0)2] 

ii) Wq = wi^2 corresponds to the fixed points w = wi^2', 

iii) For W2 < wq < wi the UV limit is w ^ W2, hence W2 is a stable fixed point of the 
theory; 

iv) For Wo > wi the UV limit is singular. The singularity is achieved at X = 1, that 
corresponds to the value 

^ 1 r/w^y/™ _ 1 

Aq O \Wo — W2'^ 
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For wq comparable to wi tg lOA,^^. We remark that the apphcabihty of the perturbative 
approach requires small Aq and therefore the singularity occurs at very high energies. For 
Wq ^ wi the position of the singularity point tg may be closer to 
zero. Since the value ts is finite, this point is singular also for ^. 
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